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1).  The  kinetics  of  two-  or  multi-component  strongly  coupled  plasmas 
requires  generalisation  of  the  methods  already  developed  by  us*  for 
one-component  systems.  This  has  been  undertaken  In  a  number  of  publications 
(A,B,C)«  Even  though  there  are  no  bound  systems  (atoms  or  Ions)  in  this  model. 


the  analysis  of  the  formalism  Is  a  necessary  prerequisite  of  the  description 

of  media  with  bound  systems.  In  (A)  we  analysed  the  cornerstones  of  our 

approximation  scheme,  the  VAA  (Velocity  Average  Approximation)  and  the  NLFDT 

(Non  Linear  Fluctuation-Dissipation  Theorem).  In  (B)  the  dynamical  NLFDT  was 

extended  to  two-component  systems.  In  (C)  we  calculated  the  nonlinear 

response  of  an  electron-ion  gas,  with  the  electron-ion  Interaction  and  quantum 

o 

effect.  Modelled  through  e  Deutsch-potentlal  . 

Publication.  (D)  and  (E)  deal  with  the  problea  of  plaaaa  oecillationa  in 
strongly  coupled  systeas.  Plata,  oecillationa  are  expected  to  reach  a  high 
level  of  excitation  (~y )  for  strong  coupling.  The  ensuing  dynaalcal 
perturbation  on  the  bound  states  bacoass  substantial  and  th«  knowledge  and 


understanding  of  the  precise  plasaoa  fttsgAMW  (mg) 
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2) .  The  standard  method  for  determining  the  degree  of  Ionization  In  a 

plasma  medium  is  through  the  use  of  the  Saha-equatlon.  However,  the 

Saha-equation,  based  on  a  model  of  a  gas  composed  of  independent  particles, 

becomes  quite  unreliable  at  high  densities*  The  main  effects  ignored  by  the 

Saha  equation  treatment  are  (i)  the  proximity  effect,  and  (11)  the  depression 

of  the  ionization  potential  due  to  screening.  We  have  developed  and  analyzed 

a  model  which  shows  the  effect  of  the  latter.  The  influence  of  the  medium  on 

the  atom  is  described  through  the  replacement  of  the  intra-atomic  potential  by 

?  1/2 

a  screened  Debye  potential  where  the  screening  constant  <  m  (4xe  n0) 
depends  on  the  free  electron  density.  Increasing  values  of  k  cause  the 
disappearance  of  high-lying  bound  states,  leading  to  the  elimination  of  the 
last  bound  state  at  ica  =  1.1  (a*  Bohr-radius)  .  Consequently,  ionization 
becomes  easier,  and  the  calculated  degree  of  ionization  is  substantially  above 
the  simple  Saha  result.  Moreover,  as  more  electrons  become  free,  the 
ionization  potential  is  further  depressed  and  an  avalanche  developes 
generating  a  sudden  phase-transltlon-llke  transition  from  the  weakly  ionized 
to  a  highly  ionized  state.  We  have  used  the  computed  tables  of  Rogers, 
Gabroske  and  Harwood^  to  describe  the  bound  states  In  a  screened  potential. 

The  Saha  equation  was  modified  by  using  the  Larkin  correction^  for  the 
partition  function.  The  resulting  rather  complex  self-consistency  condition 
for  the  degree  of  Ionization  was  solved  numerically.  Detailed  computer  work 
on  the  subject  is  in  the  process  of  completion. 

3) .  There  are  a  number  of  methods  available  for  the  study  of  the 
dynamics  of  an  atom  or  ion  Immersed  in  an  ionized  medium5.  The  Thomas-Ferml 
(TF)  method  combines  relative  simplicity  with  excellent  qualitative 
descriptive  capability.  While  the  TF  scheme  and  its  refinement  have  been 
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employed  to  the  problem  on  hand,  the  approaches  used  are  appropriate  for 
weakly  coupled  situations  and  not  for  the  case  of  strong  coupling.  The  scheme 
we  have  worked  out  is  Intended  to  Include  specific  strong  coupling  effects  and 
is  described  below. 

1.  The  system  is  considered  to  consist ,  for  the  purpose  of  the  TF 
equation,  of  three  species: 

(l)Ni  classical  nuclei  of  charge  Z,  (11)  Ne  ■  NjX  hot,  classical, 
"free”  electrons,  and  (ill)  ■  NjY  (Y-Z-X)  degenerate  bound 
electrons. 

2.  The  potential  around  a  chosen  nucleus  is  determined  by  the 
Polsson-equatlon. 


V2  Y  (x)  ■»  4re2  (Z6(x)  -  n^  (x)  -  nc(x)  +  Xn^  (x)} 
with 

nb(*)  ■  J  <*3P  1/ {  1+exp  0(€:p-u)} 
ep<0 

ep  *  P2/2m  -  Y  (x) 
ne  (*)  “  »c  U+gei^x)> 
m(x)  -  n  {1+giiOO) 

gel  aik*  Xll  ara  th*  pair  correlation  functions  for  an  electron-ion 
pair  or  ion-ion  pair. 

3.  For  the  purpose  of  the  calculation  of  the  correlation  functions  the 
system  consists  of  (1)  H*  ions  of  charge  X  and  (11)  He  hot  classical, 
"free"  electrons,  interacting  through  #M»»  *11  atk*  ♦ei* 

♦  ee  (*)  “  «2/x 


♦il  <*>  ■ 

♦•i  <*>  "  «2Z.ff,.  <*>/* 

•uch  chat  Z  >  Z«ff  >  X 
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Zeff  (o)  -  Z 

zeff  "  x 
£  being  the  ion-radius 

4,  Zeff»  i  And  Ze££,  e  are  to  be  determined  by  the  Integrated  charge 
densities  of  the  bound  electrons  only* 

5.  The  correlation  functions  gee  aiu*  8ei  are  now  to  ^ 
determined  by  using  one  of  the  strongly  coupled  static  plasma 
scheme s -either  the  STLS^  or  TI  mean  field  theories,  or  HNC. 

6*  Finally  X  is  to  be  determined  by  calculating  the  free  energy  of  the 
system  and  minimizing  it  with  respect  to  X* 

The  salient  feature  of  this  scheme  is  that  it  applies  the  philosophy  of 
the  TF-DH  method,  but  replaces  the  DH  distribution  with  the  correlation 
functions  appropriate  for  strong  coupling.  In  turn,  the  correlation  functions 
depend  on  the  effective  potentials.  A  scheme  for  the  STLS  and  TI  methods  with 
an  arbitrary  set  of  effective  potentials  for  a  two-component  system  has 
already  been  worked  out  by  us^  and  it  will  be  applied  here. 

Our  present  approach  to  the  calculation  is  to  start  with  a  constant  Zef f » 
determining  gei>  from  STLS  then  solve  for  t  and  and  then  to  Iterate  t;> 
find  a  better  Zeff*  It  will  however,  be  very  instructive  to  analyse  the 
difference  in  Y  and  nj,  between  the  values  provided  by  the  present  scheme  and 
by  those  by  the  TF-DH  scheme  as  a  function  of  y  even  for  constant  Z*ff* 
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